Nuclear symmetry energy and proton-rich reactions at intermediate energies 
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Based on an isospin dependent transport model IBUU, effects of high density behavior of nuclear 
symmetry energy in the proton-rich reaction ^'^Si+^^Si at a beam energy of 400 MeV/nucleon 
are studied. It is found that the symmetry energy affects tt^ production more than vr". More 
interestingly, comparing with neutron-rich reactions, for tt~ /tt'^ and dense matter's N/Z ratios, 
effects of symmetry energy in the proton-rich reaction both show contrary behaviors. The practical 
experiment by using the proton-rich reaction ^^Si-l-^^Ca to study nuclear symmetry energy are also 
provided. 
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Recently, pion production in heavy-ion collisions has 
attracted much attention in nuclear physics community 
[ll-Q. One important reason is that pion production is 
connected with the high density behavior of nuclear sym- 
metry energy [6*]. The latter is crucial for understand- 
ing many interesting issues in both nuclear physics and 
astrophysics 043 • The high density behavior of nu- 
clear symmetry energy, however, has been regarded as 
the most uncertain property of dense neutron-rich nu- 
clear matter [13 ■ Many microscopic and/or phe- 
nomenological many-body theories using various inter- 
actions [l^ predict that the symmetry energy in- 
creases continuously at all densities. On the other hand, 
other models 0, [l8l - [30| predict that the symmetry en- 
ergy first increases to a maximum and then may start 
decreasing at certain supra-saturation densities. Thus, 
currently the theoretical predictions on the symmetry en- 
ergy at supra-saturation densities are extremely diverse. 
To make further progress in determining the symmetry 
energy at supra-saturation densities, what is most crit- 
ically needed is some guidance from dialogues between 
experiments and transport models, which have been done 
extensively in the studies of nuclear symmetry energy at 
low densities ISll-ISal. 



Using 7r^/7r+ to probe the high density behavior of 
nuclear symmetry energy has evident advantage within 
both the A resonance model and the statistical model 
[sgI [37| . And several hadronic transport models have 
quantitatively shown that 7r"_/7r^ ratio is indeed sensitive 
to the symmetry energy [I,!!, [38, 39], especially around 
pion production threshold. These transport models, how- 
ever, usually simulate neutron-rich reactions to study the 
effect of symmetry energy, few proton-rich reactions on 
this subject were reported. Here we study the effects of 
symmetry energy on pion production, as well as the value 
of 7r^/7r+ in the proton-rich collision because the Na- 
tional Superconducting Cyclotron Laboratory at Michi- 
gan State University, Rikagaku Kenkyusho (RIKEN, The 
Institute of Physical and Chemical Research) of Japan, 
and the Cooler Storage Ring in Lanzhou, China, are plan- 
ning to do experiments of pion production to study the 
high density behavior of nuclear symmetry energy. In the 
framework of an Isospin-dependent Boltzmann-Uehling- 



Uhlenbeck (IBUU) transport model, as an example, we 
studied the effects of symmetry energy on 7r~ / 7r+ in the 
proton-rich reaction ^^Si-f-^^Si at a beam energy of 400 
MeV/nucleon. It is found that the symmetry energy af- 
fects the value of tt~ /n^ and dense matter's N/Z ratios. 
As one expected, comparing with neutron-rich reactions, 
effect of symmetry energy in the proton-rich reaction 
shows contrary behaviors. 

The isospin and momentum-dependent mean field po- 
tential used in the present work is [i^, |4l| 
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where p„ and pp denote neutron (r = 1/2) and proton 
(r — —1/2) densities, respectively. 5 = (p„ — Pp)/{pn + 
Pp) is the isospin asymmetry of nuclear medium. All pa- 
rameters in the above equation can be found in refs. [41|. 
The variable x is introduced to mimic different forms of 
the symmetry energy predicted by various many-body 
theories without changing any property of symmetric nu- 
clear matter and the value of symmetry energy at nor- 
mal density po. In this article we let the variable cc be 1. 
With these choices the symmetry energy obtained from 
the above single particle potential is consistent with the 
Hartree-Fock prediction using the original Gogny force 
(ioj and is also favored by recent studies based on FOPI 
experimental data py] . To study the effect of symmetry 
energy, we also select the stiff symmetry energy parame- 
ter X = 351 as reference. The main reaction channels 



related to pion production and absorption are 
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where R denotes A or N* resonances. In the present 
work, we use the isospin-dependent in-medium reduced 
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NN elastic scattering cross section from the scaling 
model according to nucleon effective mass [42| - |45j to 
study the effect of symmetry energy on pion production. 
Assuming in-medium NN scattering transition matrix is 
the same as that in vacuum 43:j , the elastic NN scatter- 
ing cross section in medium cr]^^™'" is reduced compared 
with their free-space value <jtf^ by a factor of 
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where ^j,nn and /i^jy reduced masses of the col- 

liding nucleon pair in free space and medium, respec- 
tively. For in-medium NN inelastic scattering cross sec- 
tion, even assuming in-medium NN — !■ NR scattering 
transition matrix is the same as that in vacuum, the den- 
sity of final states D ^ [4^ of NR is very hard to calculate 
due to the fact that the resonance's potential in matter is 
presently unknown. The in-medium NN inelastic scat- 
tering cross section is thus quite controversial (46t - i50| . 
Because the purpose of present work is just study the ef- 
fect of symmetry energy on pion production and charged 
pion ratio, to simplify the question, for the NN inelas- 
tic scattering cross section we use the free NN inelastic 
scattering cross section. The effective mass of nucleon in 
isospin asymmetric nuclear matter is 
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From the definition and Eq. ([T]), we can see that the ef- 
fective mass depends not only on density and asymmetry 
of medium but also the momentum of nucleon. 





0.10 




0.08 




0.06 








0.04 


CL 






0.02 




0.00 




affect our results here [52] ■ Please note here that, we 
do not study the effect of nucleonic distribution on the 
value of 7r^/7r+ since the values of nuclear symmetry en- 
ergy around saturation density are still open. And be- 
sides, our study here is just about the effects of high 
density symmetry energy. It is clearly shown from Fig. [1] 
that density distribution of protons is always larger than 
that of neutrons. Therefore there is no neutron skin as 
that in neutron-rich nucleus. We in the present studies 
use two proton-rich colliding nuclei 1481-1-1481 to produce 
proton-rich nuclear matter, to study the effect of high 
density behavior of nuclear symmetry energy in proton- 
rich nuclear matter. 8uch matter can form not only in the 
proton-rich nuclear reaction, but also possibly in neutron 
stars driven by the strong magnetic fields [53| . 
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FIG. 2: (Color online) Central baryon density (upper win- 
dow) and isospin asymmetry N/Z of high density region (lower 
window) for the reaction ^^Si+^^Si at a beam energy of 400 
MeV/nucleon and an impact parameter of fm. 

Fig. [2] shows the central baryon density (upper win- 
dow) and the average (N/Z)p>pg ratio (lower window) of 
all regions with baryon densities higher than pQ. It is seen 
that the maximum baryon density is about 2 times nor- 
mal nuclear matter density. Moreover, the compression 
is rather insensitive to the symmetry energy because the 
latter is relatively small compared to the E08 of symmet- 
ric matter around this density. The high density phase 
lasts for about 11 fm/c (from 2 to 13 fm/c) for this reac- 
tion. It is interesting to see that the isospin asymmetry 
of the high density region is quite sensitive to the symme- 
try energy. The soft symmetry energy (a; = 1) leads to a 
significantly lower value of [N / Z) p>p^ than the stiff one 
{x — Q). This is consistent with the well-known isospin 



FIG. 1; (Color online) Density distributions of protons and fractionation phenomenon. Because of the Esym(fi)5 



neutrons of the nucleus f|Si 

Fig. [1] shows density distributions of protons and neu- 
trons of nucleus ilSi, which were given by Skyrme- 
Hartree-Fock with 8kyrme M* force parameters [sij . 
The initializations of colliding nuclei does not evidently 
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term in the EOS of asymmetric nuclear matter, it is ener- 
getically more favorable to have a larger isospin asymme- 
try 5 in the high density region with a softer symmetry 
energy functional Esym{p)- In the supranormal density 
region, as shown in Fig. 1 of reference [54| , the symmetry 
energy changes from being soft to stiff when the param- 
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eter x varies from 1 to 0. Thus the value of {N/Z)p>pg 
becomes smaUer as the parameter x changes from 1 to 
0. It is worth mentioning that the initial value of the 
quantity {N/Z)p>pg is about 0.7 which is larger than the 
average n/p ratio of 0.57 of the reaction system. This 
is because at lower density region of proton-rich nucleus 
^|Si, there are more protons than neutrons as shown in 
Fig. [2 
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FIG. 4: (Color online) Evolution of the tt^/tt^ ratio in the 
reaction ^^Si+'^^Si at a beam energy of 400 MeV/nucleon and 
an impact parameter of fm. 



a beam energy of 400 MeV/nucleon. In the dynamics of 
pion resonance productions and decays the /i^'^)iike 
reads (mI 



FIG. 3: (Color online) Evolution of tt", 7r+ and A(1232) mul- 
tiplicities in the reaction ^^Si+^^Si at a beam energy of 400 
MeV/nucleon and an impact parameter of fm. 

To understand the dynamics of pion production and its 
dependence on the symmetry energy, we show in Fig. [3] 
the multiplicity of 7r+, 7r~ and A(1232) as a function 
of time. The multiplicity of A(1232) resonances shown 
in the figure includes all four charge states while in the 
model we do treat and follow separately different charge 
states of the A(1232). Because the value of 7r^/7r+ be- 
comes stable after 30 fm/c (can be seen from Fig. HJ, 
here we just give evolutions of A(1232) and charged pi- 
ous till 30 fm/c. At a beam energy of 400 McV/nucleon 
which is just about 100 MeV above the pion production 
threshold in nucleon-nucleon scatterings, almost all pions 
are produced through the decay of A(1232) resonances. 
The contribution due to N* resonances is negligible. It is 
interesting to see that the 7r+ multiplicity depends more 
sensitively on the symmetry energy. This is because the 
7r+ mesons are mostly produced from proton-proton col- 
lisions, where asymmetry is always larger in the reaction 
induced by the proton-rich nuclei ^^Si-|-^^Si. We can also 
see more tt^ than tt" mesons are produced. Our finding 
that in proton-rich reactions tt"*" mesons are more sen- 
sitive to the symmetry energy than 7r~ contradicts the 
results of Ref . ^] . 

To reduce the systematic errors in simulations, espe- 
cially in experimental analysis, one usually studies the 
TT instead of tt or tt^ only. Shown in 

Fig. U is effect of symmetry energy on the {t^' /iT^)iike 
as a function of time in the central reaction ^^Si-|-^^Si at 
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This ratio naturally becomes tt^/tt^ ratio at the freeze- 
out stage |5^. From Fig. |4] we can see that sensitiv- 
ity of {'iT~ /iT^)iike to the effect of symmetry energy is 
clearly shown after t — 10 fm/c. With the stiff symme- 
try (x = 0) the value of 7r^/7r+ is higher than that with 
the soft symmetry (x = 1), this is understandable within 
the statistical model for pion production [s^ . 
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FIG. 5: (Color online) The -k~ /n'^ ratio as a function of pion 
kinetic energy in the reaction ^^Si-|-^^Si at a beam energy of 
400 MeV/nucleon and an impact parameter of fm. 
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Shown in Fig. [5] is the differential 7r~/7r"'" ratios ver- 
sus the kinetic energy. In the low energy (Sfcm ^ 45 
MeV) region, around the Coulomb peak the 7r^/7r+ 
ratio is clearly separable with the x parameter varying 
from 1 to 0. Sensitivity of ■k~ /-k^ ratio to the sym- 
metry energy around the Coulomb peak is about 20%. 
In the practical experiments, the proton-rich reaction 
■^^Si-l-^^Si may be difficult to carry out. The only pur- 
pose of choosing this reaction is that we just want to show 
our studies more clearly. In the practical experimental 
plan, a reaction with proton-rich nucleus {Z > N) and 
stable nucleus {Z ^ N) (such as ^^Si-l-^°Ca) is feasible. 
Shown in Fig. |6] is the effect of symmetry energy on the 
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FIG. 6: (Color online) Evolution of the Tr~ /n^ ratio in the 
reaction ^■^Si-|-*''Ca at a beam energy of 400 MeV/nucleon 
and an impact parameter of fm. 

{tt^ /TT^)iike as a function of time in the central reac- 
tion ^^Si-|-^°Ca at a beam energy of 400 MeV/nucleon. 
Effect of nuclear symmetry energy on 7r~/7r+ is clearly 
shown. Comparing Fig. [5] with Fig. SI we can see that 
the effect of symmetry energy is smaller for ^^Si-|-^°Ca 
reaction than for ^^Si-|-^^Si. This is due to a smaller 
asymmetry (absolute) value of ^^Si-|-'*°Ca. As a com- 
parison, we also provide Fig. [71 the case of zero asym- 
metry reaction system ^^Na-|-^^Na (^^Na's half-life is 2.6 
years). As one expects, there is no symmetry energy ef- 
fect on charged pion ratio in the reaction ^^Na-t-^^Na. 



The TT /7r+ ratio is approximately 1, roughly equal to 
{5N^ + NZ)/{5Z'^ + NZ) ~ ip./Zf = (44/44)^ = 1 
in central heavy- ion reactions [5|, with N and Z being 
the total neutron and proton numbers in the participant 
region. 

In conclusion, based on an isospin dependent trans- 
port model IBUU, effects of high density behavior of 
nuclear symmetry energy on 7r~/7r^ in the proton-rich 
reactions ^^Si-|-^^Si, ^^Si-|-'"'Ca at a beam energy of 400 
MeV/nucleon are studied. It is found that the symmetry 
energy evidently affects tt"*" production and the value of 
7r^/7r+. As one expected, comparing with neutron-rich 
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FIG. 7: (Color online) Evolution of the 7r~/7r+ ratio in the 
reaction ^^Na+'^^Na at a beam energy of 400 MeV/nucleon 
and an impact parameter of fm. 



reactions, effect of symmetry energy in the proton-rich 
reaction shows contrary behavior. Studying proton-rich 
rations can not only help us to probe the symmetry en- 
ergy, but also check the theories about nuclear matter. 
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